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^ ABSTRACT 



Aims. We investigate the orf/zo-water abundance in outflows and shocks in order to improve our knowledge of shock chemistry and 



^ of the physics behind molecular outflows. 

Methods. We have used the Odin space observatory to observe the H20(lio - loi) line. We obtain strip maps and single pointings of 
^vq 13 outflows and two supernova remnants where we report detections for eight sources. We have used RADEX to compute the beam 

averaged abundances of 0-H2O relative to H2. In the case of non-detection, we derive upper limits on the abundance. 
I I Results. Observations of CO emission from the literature show that the volume density of H2 can vary to a large extent, a parameter 

that puts severe uncertainties on the derived abundances. Our analysis shows a wide range of abundances reflecting the degree to 
ry^ which shock chemistry is affecting the formation and destruction of water. We also compare our results with recent results from the 

• SWAS team. 

Conclusions. Elevated abtmdances of orf/jo-water are foimd in several sources. The abundance reaches values as high as what would 
be expected from a theoretical C-tj^e shock where all oxygen, not in the form of CO, is converted to water. However, the high 
JL abundances we derive could also be due to the low densities (derived from CO observations) that we assume. The water emission may 

^ in reaUty stem from high density regions much smaller than the Odin beam. We do not find any relationship between the abundance 

and the mass loss rate. On the other hand, there is a relation between the derived water abundance and the observed maximum outflow 
velocity. 
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1 . Introduction abundances of different species in order to distinguish between 
CO different physical scenarios. In this context, water is interesting 
Deeply embedded Class stellar systems are observed to be as- ^^^^^ that it is strongly affected by the presence of different 
0^ sociated with high velocity bipolar outflows (see e.g. Snell et al. ^p^s of shocks. At low temperatures, water is formed through a 
(3 1980) which are believed to play an important role when stars series of ion molecule reactions. This process is relatively slow 
Q>^ are formed. During the phase when material is accreted onto ^nd enhanced water abundances are thus not expected. At higher 
O the newborn star through the circumstellar disk, outflows are re- temperatures, the activation barrier for neutral-neutral reactions 
sponsible for a necessary re-distribution of angular momentum, reached, and for that reason water can be formed in a much 
. ^ The specific angular momentum of the infalhng material must more efficient way. Such elevated temperatures are reached in 
^ at some pomt decrease to allow the final collapse. Although the jow velocity shocks, where the shock is smoothed by friction 
in basic theoretical concepts can be understood, there is still a great between ions and neutrals, called continuous shocks (see e.g. 
observational need to get further knowledge about the engine of ggj-gin et al. 1998). Here, H2 is prevented from destruction and 
these flows. Different models describing the driving mechanisms enhanced water abundances are expected. In this scenario, wa- 
have been proposed and for that reason it is important to derive ter is not only formed through reactions with oxygen but can 
also be released from its frozen state on dust grains (Kaufman & 



Send offprint requests to: Neufeld 1996). In the discontinuous type of shock (jump shock), 

P. Bjerkeli, e-mail: per . b j erkeliOchalmers . se H2 is instead dissociated and water formation is prevented. The 

* Odin is a Swedish-led satellite project funded jointly by the detection of water is aggravated by the difficulty of observing 

Swedish National Space Board (SNSB), the Canadian Space Agency from ground based observatories. Prior to the launch of Odin 

(CSA), the National Technology Agency of Finland (Tekes) and Centre (j^ordh et al. 2003; Hjalmarson et al. 2003), two space born ob- 

National dEtude Spatiale (CNES). servatories capable of detecting water have been in operation, 

e n '^!?! ^'"^'^'^ Submillimetre Telescope (SEST) located at La jj^^ j^^^^^^^ ^ Observatory (ISO) (Kessler et al. 1996) and 

SiUa, Chile was funded by the Swedish Research Council (VR) and the ^, „ , ^ a ^ n 1 n-^ /otirA o\ /hi 1 • 1 

„ c .u /-vu . J ■ ■ J ■ OAAT the Submuhmeter Wave Astronomy Satellite (SWAS) (Melmck 

European Southern Observatory. It was decommissioned in 2003. 
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et al. 2000). The latter of these two also had the ability to ob- 
serve the ground state transition of 0-H2O although the beam 
size was larger (3 '3 x 4 '5 elliptical compared to 2 '1 circular for 
Odin). Water abundances derived from ISO data are in general 
higher than those derived from SWAS data. This discrepancy 
is addressed in the paper by Benedettini et al. (2002). In 2003, 
the Spitzer Space Telescope, capable of detecting hot water was 
launched (Werner et al. 2004). 

In this paper, H20(lio - loi) observations of 13 outflows 
and two supernova remnants are discussed. Shocks from su- 
pernova explosions have similar effect on the chemical condi- 
tions as molecular outflows. The different sources are discussed 
in Section (4.2) and summarized in Table 1 and 2. Table 3 in- 
cludes other outflows observed by Odin that have already been 
investigated by other authors or are in preparation for publica- 
tion (W3, Orion KL, eCha-MMSl, IRAS 16293-2422, S 140 and 
VLA1623). The analysis carried out in these papers is however 
different from the analysis made in the present paper Similar 
observations as the ones discussed here have recently been pre- 
sented by the SWAS team (Franklin et al. 2008). For that reason 
we make a brief comparison of the results for common sources. 

2. Observations and reductions 

2.1. H2O observations 

All 0-H2O observations were made with the Odin space obser- 
vatory between 2002 and 2007 (see Table 1). Each revolution 
of 96 minutes allows for 61 minutes of observations, whereas 
the source is occulted by the Earth for the remaining 35 min- 
utes. The occultations allow for frequency calibration using at- 
mospheric spectral lines. At the wavelength of the ortho-water 
ground state transition, the 1.1m Gregorian telescope has a cir- 
cular beam with Full Width Half Maximum (FWHM) of 126" 
(Frisk et al. 2003). The main beam efficiency is close to 90 % 
as measured from Jupiter mappings (Hjalmarson et al. 2003). 
The main observing mode was sky switching, where simulta- 
neous reference measurements from an unfocused 4 ?4 FWHM 
sky beam were acquired. Position switching, where the entire 
spacecraft is re-orientated in order to obtain a reference spec- 
trum, was the method of observation for a smaller number of 
targets. Three different spectrometers were used. Two of these 
are autocorrelators (ACl, AC2) and the third one is an acousto- 
optical spectrometer (AOS). The AOS has a channel spacing of 
620 kHz (0.33 kms ' at 557 GHz), while the autocorrelators can 
be used in different modes. The majority of the data have a re- 
constructed pointing offset less than 20". The data processing 
and calibration is described in detail by Olberg et al. (2003). 

3. Results 

The baseline-subtracted H2O spectra for the 15 previously 
not pubhshed sources are presented in the right column of 
Figures B.l - B.4. All spectra are smoothed to a resolution of 
0.5 km s ' . The velocity interval has been chosen in order to em- 
phasize the line profiles. In the left column, the calibrated raw 
data are plotted for comparison. These spectra have their base- 
lines subtracted using a zeroth order polynomial. 

4. Discussion 

4.1. Densities, temperatures and radiative transfer analysis 

In this paper we derive the beam averaged ortho-water abun- 
dance. The beam size is however likely to be larger than the 
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Fig. 1. The derived 0-H2O column density as a function of volume den- 
sity for different temperatures. The line intensity for this test case has 
been set to 0. 1 K while the line width has been set to 10 km s ' . 



emitting regions for several of the sources that are analyzed. In 
order to do this, we use RADEX ' (van der Tak et al. 2007). It is 
a publically available code that uses the method of mean escape 
probability for the radiative transfer For interstellar gas at rela- 
tively low temperatures and densities (<K 10^ cm"-') water exci- 
tation will be subthermal and the emission be on the linear part 
of the curve of growth, i.e. even an optically thick line will be- 
have as being effectively optically thin (Linke et al. 1977; Snell 
et al. 2000). For an effectively optically thin line, the result of a 
radiative transfer code like RADEX is essentially the same as that 
of the analytical expression for a colUsionally excited transition 

F = /iv— X„oijV(H2)yi„n(H2) „ "^"^^ (1) 

n(H2) yui 

(Liseau & Olofsson 1999). F is the integrated line flux, y\a and 
jax are the upward and downward collision coefficients, is 
the critical density and ySg is the photon escape probability. The 
molecular datafiles that are used by RADEX are taken from the 
Leiden Atomic and Molecular Database (LAMDA)^. Here, the 
collision rates between 0-H2O and H2 have been retrieved from 
Phillips et al. (1996), Dubernet & Grosjean (2002) and Faure 
et al. (2007). Recently, new rate coefficients of 0-H2O with p- 
H2 have been published by Dubernet et al. (2009). According to 
the results presented in this paper (their Figure 8), collision rates 
used in RADEX should not differ by more than a factor of three 
in the temperature regime below 40 K. The ortho to para ratio of 
H2 is in RADEX assumed to be thermal. As input parameters, the 
line intensity, line width, kinetic temperature and volume density 
of the observed gas must be supplied. The difficulty in obtaining 
good estimates for the two latter parameters has to be kept in 
mind. In Figure 1, we plot the 0-H2O column density as a func- 
tion of volume density for a test case. The line intensity has been 
set to 0.1 K and the linewidth to 10 kms From this figure, it 
is clear that the derived abundances are very uncertain when the 
kinetic temperatures are low. 



' http://www.strw.leidenuniv.nl/ moldata/radex.html 
" http://www.strw.leidenuniv.nl/ moldata/ 
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Table 1. Observation log for the sources analyzed in this paper. 



Source 


Q'(2000) 


6 (2000) 


Distance 


Backend 


Date 


'int 




(hr:min:sec) 


(deg:min:sec) 


(pc) 




(YYMMDD) 


(hr) 


L1448 


03:25:38.8 


+30:44:05 


250 


AOS 


060722 


-060813 


38 


HH211 


03:43:56.5 


+32:00:50 


315 


AOS 


070119 


- 070209 


22 


L1551 


04:31:34.1 


+ 18:08:05 


140 


AC2 


020208 


-020213 


12 


IC443-G 


06:16:43.4 


+22:32:24 


1500 


AOS 


020319 


-020321 


5 


TWHya 


11:01:51.9 


-34:42:17 


56 


AOS 


031118 


-031204 


18 


eChalN 


11:09:54.6 


-76:34 24 


150 


AC2 


031205 


-031212 


11 


Sal36 (BHR71) 


12:01:37.0 


-65:08:54 


200 


AOS 


060528 


- 060721 


41 


HH54B 


12:55:50.3 


-76:56:23 


200 


AC2 


050502 


- 050620 


10 


G327.3-0.6 


15:53:08.7 


-54:37:01 


2900 


AC2 


030731 


- 030807 


11 


NGC6334I 


17:20:53.4 


-35:47:02 


1700 


ACl 


010926 


-010927 


2 


Ser SMMl 


18:29:49.8 


+01:15:21 


310 


AC2 


041026 


-041031 


9 


3C391 BML 


18:49:22.0 


-00:57:22 


8000 


AC2 


041101 


-041106 


14 


B335 


19:37:01.0 


+07:34:11 


250 


AC2 


031024 


-031101 


7 


L1157 


20:39:06.4 


+68:02:13 


250 


AC2/AOS 


051129 


- 060127 


119 


NGC7538IRS1 


23:13:46.8 


+61:28:10 


2700 


AOS 


011214 


-011214 


1 



In this study, the temperature is taken from the literature 
while the volume density is inferred from CO observations car- 
ried out by others. These two parameters are however not ex- 
pected to have constant values across the large Odin beam due 
to the quite complex morphology of molecular outflows. For all 
the sources, we only include the cosmic microwave background 
as a radiation field. For simplicity we have chosen the line inten- 
sity to be equal to the peak value while the line width is taken 
as the width of the line at 50 % of this value. The output col- 
umn density of ortho-water is then used to derive the abundance 
relative to molecular hydrogen, X(o-H20) = A^(o-H20)/A^(H2). 



A widely recognized method to obtain the column density 
for H2 is to measure the CO abundance assuming a constant 
univers al ratio, e.g. [CO/H2] = 10""* (see e.g. Dickman 1978). 
In this paper we use this method where it is feasible (different 
methods are used for TW Hya and 3C391 BML). Volume den- 
sities and beam averaged column densities are estimated from 
literature data assuming cylindrical geometry and a mean molec- 
ular weight /i = 2.4. No correction for inclination is made. The 
inferred volume densities (Method 1) should be considered as 
lower limits for two reasons. First, the size of the water emitting 
regions is poorly known. These regions may very well be smaller 
than the CO emitting regions, potentially resulting in a higher 
average density. Secondly, shocks, if present, will compress the 
gas even further. For some of the sources, there are estimates of 
the volume density, given in the literature, that are significantly 
higher than the ones used in this paper In these cases we also 
estimate an alternative ortho-water abundance (Method 2). 



Table 2 includes the measured integrated intensity over the 
observed lines and the derived water abundance. The integrated 
intensity is measured over the entire line including the central 
region as well as the outflow wings. Excepted are those out- 
flow sources, for which strong self absorption can be seen (e.g. 
L1157, Ser SMMl). For these objects, the integrated intensity 
has been measured for the red and the blue wings separately. For 
sources with no detection, we set a 3cr upper limit on the inte- 
grated intensity in a velocity interval of 10 km s"' (except for TW 
Hya, where a linewidth of 1 km s ' has been used). 




-150:, 
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Fig. 2. The three positions observed by Odin are shown overlaid on a 
CO (2-1) map of L1448 (Bachiller et al. 1995). The circles correspond 
to the Odin beam at 557 GHz. Coordinate offsets are given with respect 
to L1448-mm: 0-2000= 03:25:38.8, (52ooo= 30:44:05.0. The positions of 
L1448-mm and L1448 1RS3 are indicated by star symbols. 

4.2. Notes on individual sources 
4.2.1. L1448 

L1448 is a dark cloud in the constellation of Perseus at the 
distance of 250 pc (Enoch et al. 2006). The large, highly col- 
limated outflow originating from L1448-mm shows enhanced 
emission from SiO in both lobes (Nisini et al. 2007). Franklin 
et al. (2008) report 0-H2O abundances of X(o-H20) = 1.5x10"^ 
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L1448 




-50 



V ["III s ] 

Fig. 3. L1448 spectra. The positions are listed in Table 2 and shown in 
Figure 2. A letter in the upper right corner indicates in which part of the 
flow the spectra have been collected (R=red, B=blue and C=center). 
The spectra have been baseline subtracted and smoothed to a resolution 
of 0.5 km s"'. 



in the blue wing and X(o-H20) - 3.7x10"*' in the red wing. We 
report observations in three positions across the structure, where 
the northern position also covers the outflow from L1448 IRS3. 
The possible detections in both lobes can, due to instabilities in 
the baselines, only be classified as likely. There is also a tenta- 
tive detection of a bullet feature in the northern position at visr ~ 
+35kms''and we note that this observation is consistent with 
the high speed CO bullet B3 reported by Bachiller et al. (1990). 
However, the preliminary analysis of HCO"^ data, recently taken 
at the Onsala Space Observatory, does not reveal any emission 
at this velocity. Mass loss rates of Mioss - 4.6x10"*' Mq yr"' for 
L1448-mm andMioss = 1.1x10"*' Mq yr"' forL1448 IRS3 were 
reported by CeccarelU et al. (1997) based on CO observations 
carried out by (Bachifler et al. 1990). A^(H2) = 6x10'^ cm"^ and 
n(H2) = IxlO"' cm"-' are inferred from mass and size estimates 
reported by the same authors. We assume the width of the flow 
to be 40". Taking the gas temperature to be T = 37 K for all po- 
sitions (Bachiller et al. 1995, dust temperature towards L1448- 
mm) we derive ortho-water abundances between 6x10""* and 
2x10"^ in the outflow. Using the higher volume density (~10'*) 
estimated by Bachiller et al. (1990), we derive ort/io- water abun- 
dances between IxlO""* and 3x10""*. 

4.2.2. HH211 

The Herbig-Haro jet HH211 is also located in Perseus, 315 pc 
away (Herbig 1998). It was observed in three dififerent positions 
enclosing the relatively small outflow. The central and northern 
beams contain the HH211-mm region. We use the mass esti- 
mates from Gueth & Guilloteau (1999) as the base for our in- 
ferred volume densities, n(H2) - 1x10^ cm~^ and column densi- 
ties, A^(H2) = 4x10'^ cm"^ in all three positions. The size of the 
outflow is taken to be 90" x 10". The temperature in the clumps 



of the IC348 region varies between 12 K close to the centers and 
20 - 30 K at file edges (BachiUer et al. 1987). Assuming a tem- 
perature T = 12 K we derive upper limits of Z(o-H20) < 8x lO""* 
at the central position and X(o-H20) < (1 - 2)xl0"^ in the out- 
flow. A temperature, T - 30K, would lower these upper limits 
by almost a factor of ten. The two-sided mass loss rate was es- 
timated by Lee et al. (2007) to be Mioss ~ (0.7 - 2.8)xl0-*' Mq 
yr~^ 

4.2.3. L1551 

LI 551 is probably one of the most rigorously studied molecu- 
lar outflows. The main source L1551 IRS5 is located at the dis- 
tance of 140pc in the Taurus- Auriga cloud complex (Kenyon 
et al. 1994). The mass loss rate is in the range 8x10"^ < Mioss 
< 2x10"^ Mq yr"' (Liseau et al. 2005, and references tiierein). 
In this paper we use the mass estimate from CO observations 
provided by Stojimirovic et al. (2006) to calculate the volume 
density and column density to n(H2) - 3x10^ cm"^ and A^(H2) 
= IxlO-^' cm"^, respectively. These authors estimate the mass of 
the outflow to be 7.2 Mq and we estimate the size of the out- 
flow to 1.3 X 0.2 pc. Assuming a kinetic temperature of 20 K 
we give upper limits on the ortho-watci abundance in three po- 
sitions across the outflow. We obtain upper limits, ranging from 
Z(o-H20) < 8x10"^ to X(o-H20) < 2x10-^ 

4.2.4. TW Hya 

At flie distance of 56 pc (Qi et al. 2008), TWHya is flie most 
nearby known T Tauri star. Its accretion disk of size 7" is seen 
essentially face-on and is point-like to Odin, i.e. particularly for 
any H2O emission region, the beam filling factor would presum- 
ably be «: 3x10"^. In our sample, the object has the lowest mass 
loss rate, viz. 10-'^'°"" < Mi„,^ <7xl0"^'' Mq yr"' (Dupree 
et al. 2005; Lamzin et al. 2004, respectively), and with an age 
of about 10 Myr (de la Reza et al. 2006), it is likely also the old- 
est. Herczeg et al. (2004) report the detection of La-excited H2, 
with an A^(H2) = 3 x lO'** cm"^, which appears accurate to within 
an order of magnitude. This value refers to the innermost disk re- 
gions. In the outer disk, CO and other trace molecules, including 
their isotopic and/or deuterated forms, have also been detected 
(Kastner et al. 1997; van ZadelhofF et al. 2001; Qi et al. 2008). 
In the HCN forming regions, the densities are determined to be 
10*' - 10** cm"^ (van Zadelhoff et al. 2001). From intermediate 
/-transitions of CO, these authors estimated temperatures to be 
in the range 40 < T < 150K. 

Not entirely unexpected, Odin did not detect the H2O 
557 GHz line^. For representative disk parameters and a line 
width of < Ikms"', the rms of 14 mK would imply an abun- 
dance, X(o-H20) < 1x10"*^. For the modehng we have used 
a temperature of 40 K. However, increasing this parameter to 
150 K wiU not decrease the derived upper limit by more than 
20%. 

4.2.5. eChaIN 

The star forming cloud Chamaeleon I is located at a distance of 
150 pc (Knude & Hog 1998). From flie estimated age (3.8xW 
yr), mass (0.2 1 Mq) and maximum velocity (~6 km s~ ' ) reported 
by Mattila et al. (1989), we obtain a mass loss rate, Mioss = 

^ Taking the difference in beam sizes into account, the Odin data are 
only a slight improvement over those obtained with SWAS (E. Bergin 
and R. Plume, private communication). 
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Fig. 4. The three positions observed by Odin are shown overlaid on a 
CO (3-2) map of Sal36 (Parise et al. 2006). The circles correspond to 
the Odin beam at 557 GHz. Coordinate offsets are given with respect 
to: Q'2ooo= 12:01:37.0, ^2000= -65:08:53.5. The positions of the sources 
IRSl and IRS2 are indicated with star symbols. 
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Fig. 5. The same as Figure 3 but for Sal 36. 

visible in Figure 4. Based on the outflow mass (1.3 Mq), dynami- 
cal time scale (1x10"* yr) and flow velocity (28 km s"') provided 
by Bourke et al. (1997)"* for the larger flows, we estimate the 
mass loss rate to be Mioss - 3.6x10"^ Mq yr"'. The wind veloc- 
ity is assumed to be 100 kms"'. 



3.3x10 ^ Mq yr The velocity of the wind is assumed to be 
100 kms-i. We estimate NiHi) = 4xl02"cm-2 and n(H2) = 
2x10"' cm""* from CO observations carried out by the same au- 
thors. The width of the flow is approximately 0.1 pc. Adopting 
the temperature 50 K, given by Henning et al. (1993), we obtain 
an upper limit, ^(o-HaO) < 3x10"^. 

4.2.6. Sa136 (BHR71) 

The Sal36 Bok globule outflows (Sandqvist 1977), at about 
200 pc distance, are driven by a binary protostellar system. The 
secondary CO outflow driven by IRS 2 is more compact (Bourke 
2001) than the larger outflow driven by IRS 1. 

Taking the column density measurements made by Parise 
et al. (2006), assuming a [CO/H2] ratio of lO""* we obtain A^(H2) 
= IxlO^i cm-2, N(U2) = 4x10^" cm^^ and A?(H2) = 3x10^1 cm^^ 
in the red, central and blue part of the flows respectively. We as- 
sume that the flow has a depth of 0.07 pc, yielding volume den- 
sities n(H2) = 6x10"^ cm-\ n(H2) = 2x10^ cm^^ and n(H2) = 
IxlO'^cm"-' in the same regions. Parise et al. (2008) estimate 
r = 30 - 50 K from CO and methanol observations. In our mod- 
eling we use r = 40 K. We estimate the ortho-water abundances 
to (0.1 - l)xlO"^ in the outflow and 2x10""* at the central po- 
sition. However, Parise et al. (2008) give a density of «(H2) - 
1x10^ cm"^ in the region. Using this higher value we obtain or- 
tho-water abundances of (2 - 6)xl0"^ in the outflow and 3x10"* 
towards the central source. The emission has broader wings in 
the central position, a feature present also in the SWAS data. The 
origin of this high velocity component and the elevated water 
abundance might be the smaller outflow originating from IRS 2, 



4.2.7. HH54B 

The Herbig-Haro object HH54B is situated in the Chall cloud 
at roughly 200 pc (Hughes & Haitigan 1992). During the obser- 
vations, the correlator suffered from ripple. This had the effect 
of increasing the line intensity and a substantial amount of data 
had to be abandoned. The signal detected can therefore only be 
classified as tentative although we are confident that the data do 
not show any systematic variations. Complementary data were 
obtained in CO(3-2), CO(2-l), SiO(5-4), SiO(3-2) and SiO(2- 
1) with the SEST telescope (see Appendix A) and in CO(5-4) 
with Odin. The CO(5-4) data suffer from frequency drift, some- 
thing that puts an ~10kms"' uncertainty to our velocity scale. 
However, we do not believe that this gives large uncertainties on 
the line strength. 

No shock-enhanced emission was detected in any of the ob- 
served SiO transitions. This result seems not easily reconcilable 
with the prediction from theoretical C-shock models (see Fig. 6 
of Gusdoif et al. 2008), which appear closely adaptable to the 
conditions in HH54 (Liseau et al. 1996; Neufeld et al. 2006). 

Both the CO (2-1) and (3-2) lines show an absorption fea- 
ture at -1-2.4 km s"', which corresponds to the LSR-velocity of 
the molecular cloud. In addition, a strong blue wing is observed 
in all positions, but essentially no redshifted gas, which is in 
agreement with the (1-0) observations by Knee (1992). Both (2- 
1) and (3-2) transitions peak at the central map position, i.e. on 
HH54B itself and their integrated intensities, Jt* dv, are given 
in Table A.l. 

These authors estimate to 10 700 years and 10 200 years for the 
south east and the north west lobe respectively. 



6 



P. Bjerkeli et al.: Odin observations of water in molecular outflows and shocks' 



HH54 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 

C0(5-4) 








, 1 , , 1 1 \ , 1 , , \ K 

C0(3-2) A 


11 1 11 , 




— -^^■"^"-^^ 


— , — 1 — , — , — 1 — 1 — \ — , — 1 — , — , — \ — h 

; co(2-i) 1^ 


— \ — 1 — 1 — — \ — 1 — \ — : 

A : 
-^^^^ : 


1 .... 1 












1 . . , , 1 . , . . 1 . 





-20 -10 



Vi„ [to s-l] 

Fig. 6. From top to bottom the CO(5-4), CO(3-2), CO(2-l) and 
H20(lio - Iqi) spectra observed towards HH54B are plotted. The dot- 
ted vertical line shows the cloud LSR- velocity at +2.4 km s"' . 

For the comparison with the ISO-LWS model of Liseau et al. 
(1996), we use the average radiation temperature, approximated 
by < 7- >= 77-^ fr* dvl Jdv. This yields < T21 >= (1 .2 ± 0.02) K 
and < >= (3.2 + 0.33) K for the CO (2-1) and (3-2) lines, 
respectively. Both values are smaller, by 25% and 40% respec- 
tively, than the model predictions of 1 .6 K and 5.7 K, which were 
based on a single- temperature approximation. Even though the 
strength of the CO(5-4) line is uncertain, it shows a slightly 
higher temperature than the predicted 0.3 K. Putting it all to- 
gether, the model predicts the radiation temperature in all three 
lines to within an order of 2. For the RADEX analysis we use the 
ISO-LWS model T = 330 K, n(H2) = 2xl05cm-3 and NiUi) 
= 3x10^' cm"-^, where the column density has been diluted to 
the Odin beam. We obtain a beam averaged water abundance of 
Xio-HiO) = 3x10 ^ However, recently Neufeld et al. (2006) 
estimate a higher H2 column density for the warm gas, a fact 
that could alter our derived abundance by an order of magnitude 
downwards. The mass loss rate of the unknown driving source 
has been estimated by Giannini et al. (2006) as Mioss = 3 X 10"^ 
Mq yr-i. 

4.2.8. G327.3-0.6 

The hot core G327.3-0.6 is located in the southern hemisphere 
at the distance of 2.9 kpc (Bergman 1992). CO line profiles ob- 
tained by Wyrowski et al. (2006) weakly indicate the presence of 
outflows, however, to date there is no further study of this. From 
CO observations performed by these authors, we infer A'(H2) 



NGC6334I 
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Fig. 7. The same as Figure 3 but for NGC6334 1. The black spectrum 
represents the first case and the grey spectrum represents the second 
case as described in Section 4.2.9 



= 2x10^2 cm-2 and n(H2) = 4x10^ cm-^ assuming ["^CO/H2] 
= 10"^ and a source size 25". The volume density inferred is 
slightly lower than the range 10^-10^ cm~^ given by Bergman 
(1992), who also estimates the temperature to be within the 
range of 40 - 200 K. Using a value of 100 K gives an upper limit 
of X(o-U20) < 8x10"'°. A temperature of 40 K will increase the 
inferred upper limit by a factor of 8. 

4.2.9. NGC6334I 

At least two outflows are emerging from NGC6334 1, located in 
the constellation Scorpius (McCutcheon et al. 2000) at the dis- 
tance of 1.7 kpc (Neckel 1978). From CO observations provided 
by Leurini et al. (2006) we obtain a beam averaged column den- 
sity 1x10^" cm"^ and a volume density 4x10^ cm"^. We assume 
that the gas temperature is the same as the dust temperature, viz. 
T = 100 K (Sandell 2000). This is consistent with Leurini et al. 
(2006) who set a lower hmit on the kinetic temperature at 50 K. 
With the above properties we derive an abundance of X(o-H20) 
= 5x10"^ 

The baseline subtracted spectrum does not show any evi- 
dence of high velocity gas. However, we do not find this easily 
reconcilable with the high velocity gas detected in several CO 
transitions by Leurini et al. (2006). One possibiUty could be the 
curved baseline hiding the outflow wings. Therefore, we inves- 
tigate also an alternative case where we assume that the entire 
curvature stems from the outflowing gas. This secondary sce- 
nario is perhaps not very likely. However, at present it is not 
possible to draw any firm conclusions. The estimated depth of 
the absorption feature is greater than the continuum level of 360 
Jy, interpolated from 800 yum observations provided by Sandell 
(2000). In this secondary case we derive the abundance X(o- 
H2O) = 2x10-^ 

4.2.10. SerSMMI 

The Serpens star forming dark cloud is situated in the inner 
Galaxy at a distance of 310pc (de Lara et al. 1991). Franklin 
et al. (2008) estimated ortho-water abundances of X(o-H20) = 
7.1x10"^ and ^(o-HzO) = 3.8x10"^ in tiie blue and red wing 
respectively while Larsson et al. (2002) estimate the water abun- 
dance to X(U20) - 1x10"^ in the region. Davis et al. (1999) 
provide the mass and size of the outflow based on CO observa- 
tions. Assuming a width of 0.2 pc yields n(H2) = 1x10^ cm~^ 
and N{H2) = 5x10^*^ cm~^. From Ol(63/im) measurements car- 
ried out by Larsson et al. (2002, and references therein) we ob- 
tain the mass loss rate, Mjoss = 3x10"^ Mq yr"'. The temperature 
of the dust was constrained by White et al. (1995) to be 30 K 
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Fig. 8. The same as Figure 3 but for Ser SMMl. 



< r < 40 K. Using T — 35K and the above properties we obtain 
X(o-H20) = 9x10-5 andZ(o-H20) = 5x10 ^ in the blue and red 
flow. 

4.2.11. B335 

The dense core in the B335 globule is believed to be one of the 
major candidates for protostellar collapse. This isolated source at 
the distance of 150pc (Stutz et al. 2008) harbors several Herbig- 
Haro objects associated with a bipolar outflow (see e.g. Galfalk 
& Olofsson 2007). Following the outflow mass estimate of 0.44 
Mq in Hirano et al. (1988) we derive n(H2) = IxlO"' cm"^ and 
A^(H2) = 3x10^° cnT^. The size of the outflow is taken to be 2' 
X 8'. Based on the same mass, a dynamical time scale (2.3x10"* 
yr) and a flow velocity (13 km s"'), provided by these authors we 
derive Mioss - 2.5x10"'' Mq. The velocity of the wind is assumed 
to be 100 kms"'. In our modeling we use a kinetic temperature 
of 20 K. This is an intermediate value in the modeling carried 
out by Evans et al. (2005, their Figure 8) who give the temper- 
ature as a function of radius for the inner 0. 1 pc region. The de- 
rived upper limit of the water abundance is, X(o-H20) < 1x10"^. 
The rms implies a water abundance, Xio-HjO) < ~ 10"^ accord- 
ing to the infall model presented by Hartstein & Liseau (1998). 
Ashby et al. (2000) estimate the upper limit on the ortho-water 
abundance to be within the range 1.3x10"^ - 7x10"^. 

4.2.12. L1 157 

LI 157 is a Class object in the constellation of Cepheus sit- 
uated at the distance of 250 pc (Looney et al. 2007). It drives 
a prototype bipolar outflow which was observed by Odin in 
four different positions. The mass loss rate can be derived from 
mass (0.62 Mq) and time scale (15 000 yr) estimations given 
by Bachiller et al. (2001). Assuming a maximum CO velocity 
of 20 kms"' and a stellar wind velocity of 100 km s"', we de- 
rive Mioss - 8.3x10"^ Mq yr"' . The Odin strip map covers the 
bulk of the outflow with two pointings in the red wing, one in 
the blue, and one on the driving source of the outflow itself 
(Figure 9). The observations were carried out using the AOS 
and the AC2 simultaneously. The spectra shown in Figure 10 
are the averages of the merged data. The total mass in the dif- 
ferent parts of the flow were obtained from CO observations 
carried out by Bachiller et al. (2001). From this we estimate 
n(H2) = 2x10^ cm"^ and A^(H2) = 2x102" cm"^ in the northern 
lobe, while n(H2) = 3xl03cm-'' and N(H2) = 2x102° cm^^ in 
the central and southern region. The size of the CO outflow is 
taken to be 50" x 375". For all four positions we set the ki- 
netic temperature to T - 30 K, a rough global estimate based 
on Bachiller et al. (2001). We calculate the abundances in the 
outflow to be within the range of 2 x lO"'* and 1 x 10"^. The 




100 50 -50 -100 
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Fig. 9. The four positions observed by Odin are shown overlaid on a 
CO (2-1) map of LI 157 (Bachiller et al. 2001). The circles correspond 
to the Odin beam at 557 GHz. Coordinate offsets are given with re- 
spect to L1157-mm, indicated in the figure with a star symbol: 0-2000= 
20:39:06.4, ^2000= -1-68:02:13.0. The black and white squares refer to 
the knots BO, Bl, B2, RO, Rl, R and R2 described in Bachiller et al. 
(2001). 
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Fig. 10. The same as Figure 3 but for LI 157. 
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derived water abundance in the central region is slightly lower, 
X(o-H20) = 2x10 in the blue lobe and ^(o-HaO) = 3x10"^ 
in the red. The increased blue emission Ukely originates from 
the outflow. Within the Odin beam are the positions BO and Bl 
that show peaked emission in H2CO, CS, CH3OH, SO (Bachiller 
et al. 2001, their Figure 1) and SiO (Nisini et al. 2007). Frankhn 
et al. (2008) estimated XCo-HjO) = S.OxlO""^ and XCo-HjO) = 
9.7x10"^ in the blue and red wing respectively. 

Bachiller et al. (2001) estimates the density around the pro- 
tostar to be ~10''cm"^. When moving from BO to B2, the den- 
sity changes from ~3 to 6x10^ cm"-'. Using n{H2) - 1x10*" cm"^ 
and n(H2) = 5x10^ cm"^ for the central and southern part re- 
spectively we obtain the abundances ^(o-HiO) = 5x10"^ and 
X(o-H20) = 2x10"^. 

4.2.13. NGC7538 IRSl 

NGC7538 IRSl is a region of ongoing high mass star formation. 
The main infrared source IRS 1 is located at the boundary of an 
HII region in the Perseus arm, located at a distance of 2.7 kpc 
(MoscadeUi et al. 2008). In addition to IRSl, and its high ve- 
locity outflow, also several other sub-mm sources fall into the 
large Odin beam. The mass loss rate from IRS 1 was estimated 
by Kameya et al. (1989) to be Mioss = Ix 10 Mq yr ' . The total 
mass and size of this region is given by the same authors, yield- 
ing a volume density and column density of n(H2) = 3x10^ cm"^ 
and A^(H2) = 4x10^*^ cm~^ respectively. The kinetic temperature 
is assumed to be the same as the dust temperature, i.e. T = 40 
(Sandell & Sievers 2004). We infer a Utnit to the water abun- 
dance of X(o-H20) < 1x10-^ 

4.2.14. Supernova remnants 

3C391 BML is a supernova remnant at the distance of 8 kpc 
(Chen & Slane 2001). The temperature, volume density and col- 
umn density of the gas were constrained to 50 K < T < 125 K, 
10"^ cm-3 < n(H2) < 5x10^ cm ^ and N(U2) > 4x10^0 cm~^ from 
OH 1720 MHz maser observations carried out by Lockett et al. 
(1999). Using the lowest values for the temperature and volume 
density we calculate an upper limit of X(o-H20) < 4x10"^. 
IC443 is a shell type supernova remnant at the distance of about 
1.5 kpc (Fesen 1984). The broad emission peaks are labeled A 
through H, where Odin has observed clump G. The volume den- 
sity and temperature were modeled by van Dishoeck et al. (1993) 
as n(H2) = 5x10^ cm^^ and T = 100 K. Using the CO column 
density inferred by the same authors, assuming a [CO/H2] ratio 
of lO""* and a source size of 40" x 100", we obtain N(ll2) = 
3x10^^ cm-2. The derived abundance is X(o-H20) = 4 x 10"^ 
This is in agreement with Snell et al. (2005) who derive an o- 
H2O abundance with respect to '^CO, X(o-H20) = 3.7x10-"*. 
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Fig. 12. The same as Figure 3 but for IC443-G. 

4.3. Water abundance 

The water abundances inferred from our analysis are given with 
respect to the molecular column density within the Odin beam. 
The values span a wide range with most of the sources having 
an abundance of the order lO^^-lO""*. The highest abundances 
in our sample are found in the outflows of LI 157 and L1448. 
Supposing that the assumed physical properties are correct, we 
see an increased abundance in both the blue and the red wing of 
LI 157. The derived water abundances are significantly low in 
G327. 3-0.6. The reason for this is the combination of the high 
H2 column density and temperature used. The water emission is 
also very Ukely beam diluted due to the large distance, resulting 
in lower beam averaged 0-H2O abundances. 

Assuming an inclination angle of 60° with respect to the 
line of sight for all the targets we plot the 0-H2O abundances 
versus the maximum velocities (Figure 13). The maximum 
velocities are taken from the spectra as the maximum offsets 
between the cloud velocity and the flow velocity. There is a 
correlation between the derived abundances and the maximum 
velocities of the outflowing gas. The solid line in Figure 13 is 
the first order polynomial least square fit: 

Z(o-H20) = lO-'vl,,. (2) 

The correlation coefficient is 0.57 while the p- value, testing the 
hypothesis of no correlation is 0.04. Following the C-shock mod- 
elling carried out by Kaufman & Neufeld (1996), a relation- 
ship is expected when shocks are responsible for the emission. 
However, the abundances do not show any tendency to level off 
with velocities larger than ~20 km s"' as shown in their Figure 3. 
This is an indication of a lack of /-type dissociative shocks 
in this sample of outflows. There are however uncertainties in 
our simplified analysis, for example the inclination of the out- 
flow lobes and the RADEX analysis. The dependence between 
the water abundance and the maximum velocity is consistent 
with the analysis made by Franklin et al. (2008), although their 
method is based on binning the weak outflow lines into inter- 
vals of 5 km s"' . We also plot the ortho-water abundances ver- 
sus mass loss rates in Figure 14. G.327-0.6 is not included due 
to the lack of observations towards the assumed outflow, neither 
is NGC6334I included in Figure 13 - 14 due to the curvature in 
the baseUne. No obvious relationship can be picked out. The ab- 
sence of a correlation is surprising inasmuch as high mass loss 
rates affect the budget of material available for water formation. 

4.4. Comparison with SWAS data 

The gas volume density, used as an input parameter to RADEX , 
has for all the outflows except TW Hya been derived from CO 
observations. This method generally underestimates the mass of 
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Fig. 13. 0-H2O abundance plotted against the maximum velocity for an 
overall inclination of 60° with respect to the line of sight. The solid 
line represents a linear fit of log[X(c)-H20)] versus log(i'm!a) with the 
same inclination angle applied. The dotted dashed line represents the fit 
with a inclination correction of 35° while the dashed Une represents the 
fit with an inclination angle correction of 85°. The errorbars from the 
measurements are smaller than the circles. 
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the regions. The possibility that the water emission is originat- 
ing from gas with a higher voltime density than this can therefore 
not be ruled out. Nevertheless, we are confident that the volume 
density is actually spanning over a wide range of values. This is 
also one of the reasons why several of our derived abundances 
deviate from those inferred by Franklin et al. (2008). They use a 
single volume density of niWj) = 10^ cm"^ and our values differ 
from this by more than two orders of magnitude for some of the 
sources (see Table 2). Figure 15 shows a histogram of the num- 
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Fig. 15. A histogram of the gas volume density estimated in the outflows 
studied in this paper shows a variation that spans over six orders of 
magnitude. 



bers of sources within different volume density ranges. The dif- 
ficulty to estimate the gas column density of the water emitting 
regions is a problem that has to be adressed in order to interpret 
the future observations with Herschel. 

The integrated intensities for the SWAS and Odin outflow 
spectra are compared in Figure 16 to provide an estimate of the 
source size of the water emitting regions. The dashed 1:1 ratio 
line illustrates the case where the source is filhng both antenna 
beams while the dashed 3.37:1 line is indicating a small source 
size compared to both beams. Assuming that both the emitting 
sources and the antenna responses are circularly symmetric and 
Gaussian, these ratios follow from the relation: 



/odin _ (3 '3 X 60 X 4^5 X 60) -1- 6'^ 
hwAs 



1262 -I- 6*2 



(3) 



where 6 is source size. The SWAS spectra have been retrieved 
from the SWAS spectrum server and baseline subtracted. The 
integrated intensities are measured in the same regions as for the 
Odin spectra and the spectra chosen for comparison are those 
where the pointing is equal or within a smaller fraction of a 
beam. We can conclude that the sources are filling a large frac- 
tion of the beam for Ser SMMl, LI 157 and NGC 7538 1 . On the 
other hand, the source sizes are small for HH54 B and L1448. 

4.5. Outflows and observed water abundances 

The main objective for these observations was to use the ground 
state 0-H2O transition as a tracer for shocked gas. Available 
shock models by Bergin et al. (1998) show that a shock velocity 
in excess of lOkms"' should convert essentially all oxygen, not 
in the form of CO, into water resulting in water abundances rela- 
tive to H2 of the order 10 This water is also believed to persist 
for ~10^ years. For several of our sources, the water abundance 
increases as high as what would be theoretically possible for 
these C-type shocks. The fact that we are observing abundances 
as high as this, might be an indication that the shock velocities 
are in this regime. On the other hand, the large Odin beam does 
not show whether the water emission regions are small or not. 
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Fig. 16. Integrated intensities for common sources of Odin and SWAS. 
The dashed lines show the 3.37:1 and the 1:1 ratios between the Odin 
and SWAS integrated intensities. The error bars refer to the analysis and 
the solid lines represents ratios 1:1 and 3.37:1 with a 15 % uncertainty 
applied. This is the estimated error limit from the data reduction. 

The density in these regions may also very well be orders of 
magnitude higher than the values used in this paper, a fact that 
could alter our inferred values significantly. Future observations 
with the Herschel space observatory will give us the ability to 
spatially resolve outflows and it will also give us the possibil- 
ity to observe several different transitions originating in gas of 
different temperatures. This will hopefully give a much deeper 
insight into the chemistry occurring in these structures. 

5. Conclusions 

We make the following primary conclusions: 

1. We have observed 13 outflows and two supernova remnants 
and detect the ortho-water ground state rotational transition 
in seven outflows and one super-nova remnant. 

2. The column densities of 0-H2O have been investigated with 
RADEX , having the volume density, temperature, line inten- 
sity and linewidth as input. Elevated abundances of water are 
found in several sources. The abundances are as high as one 
would expect if all gaseous oxygen had been converted to 
water in a C-type shock. 

3. There is no distinct relationship between the water abun- 
dance and the mass loss rate. 

4. There is a correlation between the 0-H2O abundance and the 
maximum velocity of the gas. 
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Table 3. Column densities of 0-H2O and estimates of the ortho-v/ater abundance, X(o-H20) = iV(o-H20)/iV(H2). 



Source 


Offset / T^(H20)dv 
(",") (Kkms-i) 


All 


(iU cm ) 


(lU cm ) 


A(Cl-hl2U) 

(10 ) 


Odin sources analyzed by other authors: 










W3« 


(0,0) 3.9(0.1) 








0.02 


Orion KL* 


(0, 0) 323 






400 


80 / 1000 




(+0, -240) 24 






0.09 


0.1 


Cha-MMSP 


(+18, -21) 








<0.07 




(-35, +32) 








<0.07 


IRAS 16293-2422'' 


(+68, -50) 


red wing 






30 






blue wing 






11 




(+31, -10) 


red wing 






600 






blue wing 






60 




(-23, -8) 


red wing 






130 






blue wing 






11 


8140" 


(0, 0) 








0.05 - 4 


VLA1623^ 


(0, 0) 


red wing 


2 


0.2-0.4 


100 - 200 



Notes to the Table: " Wilson et al. (2003). The offset is with respect to 0-2000 = 02:25:40.6, ^2000 = +62:05:57. * Olofsson et al. (2003). These 
authors use two different values for Af(H2) towards the (0, 0) position. The offsets are with respect to a2ooo = 05:35:14.4, d^im = -05:22:30. Klotz 
et al. (2008). The offsets are with respect to Qr2ooo = 11:06:31.7, (^2000 = -77:23:32. RistorceUi et al. (2005); RistorceUi, private communication. 
The offsets are with respect to ^2000 = 16:32:23.0, ^2000 = -24:28:40. ' Persson et al. (2009), towards the PDR and the outflow. The offset is with 
respect to a2ooo = 22:19:19.4, ^2000 = +63:18:50. ^ Larsson et al. in prep, private communication. The offset is with respect to a2ooo = 16:23:20.2, 
.52000 = -24:22:55. 



Appendix A: Ground based observations 

The observations with the Swedish ESO Submillimetre 
Telescope (SEST) were made during 11-20 August 1997 and 2- 
6 August 1998. Some complementary SiO(2-l) map data were 
collected during 7-9 February 2003. The observed molecules 
and their transitions are listed in Table A.l and for SiO (2-1) and 
(3-2), the observations were done simultaneously. 

SIS receivers were used as frontends and the backend was 
a 2 X 1 GHz multi-channel acousto-optical spectrometer (AOS), 
split into two halves for the SiO (2-1) and (3-2) observations. 
For the SiO (5-4) and the CO observations, both the high reso- 
lution (HRS) and the low resolution (LRS) backends were used. 
The channel separation of the HRS was 43 kHz and the spec- 
tral resolution 80 kHz covering the bandwidth of 86 MHz, cor- 
responding to 6v = 0.07 km s"^ and Av = +35kms ' for the 
velocity resolution and range at 345 GHz, respectively. For the 
LRS, these parameters were channel separation 0.69 MHz and 
Av = 1.4 MHz spectral resolution, providing 6v = 1.2kms~' 
and Al" = ±700kms~\ respectively, for the CO (3-2) observa- 
tions. 

The data were chopper-wheel calibrated in the r*-scale 
(Ulich & Haas 1976) and the main beam efficiencies at the dif- 
ferent frequencies, rj^^^,, are given in Table A.l. The pointing of 
the telescope was regularly checked towards point sources, mas- 
ing in the SiO (v=l, J =2-1) line, and was determined to be better 
than 3" (rms). However, for HH54, all SiO data refer to the vi- 
brational ground state, v=0, and the (2-1) and (3-2) data were 
obtained in frequency switching mode, with a frequency chop of 
4 MHz. 

Knee (1992) assigned a kinetic gas temperature of ~ 15 K to 
the bulk cloud material. In the CO (3-2) line, this should yield 
a high contrast between the cloud and the high velocity gas. To 
achieve flat, optimum baselines, the CO (3-2) observations were 
done therefore in wide dual beam switching with throws of ±1 1' 
in azimuth. At this maximum amplitude available at the SEST, 
the reference beams were still inside the molecular cloud. For 



this reason, this mode could not be adopted for the CO (2-1) 
observations, which were performed in total power mode. The 
reference position was 1° north of HH54. Centered on the ob- 
ject, nine point maps with 25" spacings were obtained in both 
CO lines. In addition, a tighter sampled nine point map with 15" 
spacings was also made in CO (3-2). 
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Table A.l. Molecular line observations with the 15 m SEST 



Transition 


Frequency 


HPBW 


Imb 


^sys 


^rms 


Backend 








(MHz) 


(") 




(K) 


(mK) 




(kms-') 


(Kkms-i) 


CO (2-1) 


230538.0 


23 


0.50 


500 


95/10 


HRS/LRS 


-26.3 to -0.3 


15.6 ±0.3 


CO (3-2) 


345796.0 


15 


0.25 


1600 


470/70 


HRS/LRS 


-16.4to-F0.8 


13.6 ± 1.4 


SiO (2-1) 


86846.9 


58 


0.75 


160 


10" 


HRS 






SiO (3-2) 


130268.7 


38 


0.68 


190 


15" 


HRS 






SiO (5-4) 


217104.9 


23 


0.61 


750 


35/7" 


HRS/LRS 







Note to the Table: " This is the rms of the average spectrum, including all positions. 
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Appendix B: Online material 

For comparision we show the calibrated rawdata together with 
the baseline subtracted and smoothed spectra already shown in 
the text. 
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Fig.B.l. This figure shows the L1448, HH211 and 1C443-G spectra. The positions are listed in Table 1. All spectra on the right are smoothed to 
a resolution of 0.5 kms"' and clipped for clarity. These spectra are also calibrated in frequency. The spectra on the left are the calibrated rawdata 
with a zero order baseline subtraction. For strip maps, a letter in the upper right comer indicates in which part of the flow the spectra has been 
collected (R=red, B=blue and C=center). 
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Fig.B.2. The same as Figure B.l but for L1551, Sal36 and TW Hya. 
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Fig.B.3. The same as Figure B.l but for eps ChaIN, HH54 B, G327.3-0.6, NGC6334I and Ser SMMl. 
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rig.B.4. The same as Figure B.l but for 3C391 BML, B335, NGC7538 IRSl and L1157. For L1157, the AC2 data are plotted in gray and the 
AOS data in black. 



